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We examine the influence of coupUng between particle-hole and particle-particle spin fluctuations 
on the inelastic neutron scattering resonance (INSR) in cuprate superconductors in both weak and 
strong interaction limits. For weak-interactions in the particle-hole channel, we flnd that the inter- 
channel coupling can eliminate the resonance. For strong interactions which drive the system close 
to a Q = (tt, tt) magnetic instablity, the resonance frequency always approaches zero but its value 
is influenced by the interchannel coupling. We comment on constraints imposed on cuprate physics 
by the INSR phenomenology, and a comparison between the cuprates and the newly-discovered iron 
pnictide superconductors is discussed. 

PACS numbers: 



I. INTRODUCTION 

The low-temperature properties of cuprate supercon- 
ductors seem to be well described by mean-field-theory 
with effective interactions which lead to d-wave super- 
conductivity, except possibly in the extreme underdoped 
limit. At the same time the apparently ubiquitousi pres- 
ence of a low- frequency spin- resonance near Q = (tt, tt) 
in inelastic neutron scattering experiments and the ap- 
pearance of spin order in some systems^ when an exter- 
nal magnetic field is applied, suggest that superconduct- 
ing cuprates are close to an antiferromagnetic instabil- 
ity. In previous work2, we have argued that coupling^ be- 
tween particle-hole and particle-particle spin fluctuations 
likely contributes to the suppression of the superfluid 
density^ observed in underdoped cuprates. In this arti- 
cle we specifically address, within time-dependent mean- 
field theory, the influence of particle-particle to particle- 
hole coupling on the inelastic neutron scattering reso- 
nance (INSR) position. When the interchannel coupling 
is ignored^ an arbitrarily weak particle-hole channel in- 
teraction will induce a resonance. We flnd that the in- 
terchannel coupling can eliminate weak-interaction res- 
onances. For stronger interactions the system can be 
driven close to an antiferromagnetic instability. In this 
limit we find that the way in which the resonance fre- 
quency approaches zero is influenced by the interchannel 
coupling. There is no limit in which RPA-type theories 
without the interchannel coupling yield reliable resonance 
frequency predictions. 



II. RELATIONSHIP BETWEEN THE 
DIAGONAL AND OFF-DIAGONAL RPA 
SUSCEPTIBILITIES 

In time-dependent mean-fleld theory the most general 
linear response of a system consists of correlated particle- 
hole and particle-particle excitations. For square-lattice 
models with on-site and near-neighbor electron-electron 
interactions we have shown previously^ that at Q = 



(tt, tt) correlations occur between between four different 
weighted particle-hole transition sums. The most general 
response function can therefore be evaluated by invert- 
ing four by four matrices. For a model (motivated by the 
cuprate superconductivity literature) with only on-site 
([/) and near-neighbor [V) spin-independent interactions 
and nearest-neighbor Heisenberg (J) spin interactions we 
flndZ that: 
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where V = diag(-C/ - 2 J, J/2 -2V,V + J/4, V + J/4) is 
the interaction kernel. 
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is the bare mean-fleld-quasiparticle response function, 
and k' = k + Q reduced to the flrst Brillouin-zone. The 
indices a and b refers to the four coupled two-particle ex- 
citation channels and fa is the corresponding coherence 
factor deflned in Ref. Q- As shown in previous worki 
the INSR mode is dominated by two channels, the spin 
flip {a — I) and the c?-wave pair phase (a = 4) modes, 
when parameters are in a range broadly consistent with 
experiment. The same conclusion can be reached on the 
basis of earlier approximate RPA theories^ which omitted 
the second channel. The qualitative discussion of INSR 
physics below is based on a model in which only these two 
modes are retained, although the full 4-channel model is 
required for quantitative accuracy. 

In the truncated two-channel model the INSR solves 



= KM - Vs 



C^{uj)lu^ (3) 



where K'^iuj) 

Vs = U+2J, Vcf, — V + J/A, and the frequency-dependent 
stiffness and Berry-phase coupling parameters are defined 
by: 
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K^{bj = 0) and K'^(uj — 0) are"^ proportional to the 
energetic cost of static spin-density and spin-polarized 
electron-pair density fluctuations, as estimated by mean- 
field-theory while C{uj = 0) is^ the Berry-curvature asso- 
ciated with adiabatic evolution of electron-electron and 
electron-hole pair mean-fields. Stability of the d-wave su- 
perconducting state requires that both K^{uj = 0) and 
K'^{uj = 0) be positive. In the two-channel approxima- 
tion the interchannel coupling is captured by the function 
Xi4.9p('^)- If Xi4,(ip(<^) was negligible, Eq. [3] would be sat- 
isfied at both the spin-exciton resonance frequency tOex 
defined by 

1/Xll,gp(^e2;) +Vs = 0, (5) 

and at the 7r-resonance frequency cj^ defined by 

1/X44,<}p(w7r) - = 0. (6) 

When coupling between particle-hole and particle- 
particle channels is neglected both -R(w) and the Berry 
phase coupling C((jj)'^'^ are set to zero. It is therefore 
important to study how R{uj) and C(w) depend on the 
model parameters. Since both R and C are dimension- 
less, the unbiased justification of whether or not the in- 
terchannel coupling is negligible should be based on these 
two quantities. It is evident from the definition of R{uj) 
and from the way in which it appears in the above equa- 
tions, that the interchannel coupling certainly can not 
be neglected when it approaches 1. To judge the im- 
portance of the Berry-phase coupling C{u!), we should 
compare |C| x cOres with typical values for ^yK^xK^, 
which indicates that values of |C| comparable to 1 also 
correspond to strong interchannel coupling. 

The mean-field-theory of INSR is defined by the hop- 
ping parameters of the square-lattice tight-binding model 
and by the interaction parameters U,V, and J. We em- 
phasize that all these parameters should be viewed to a 
large degree as effective parameters which depend on mi- 
croscopic many-body physics in ways which are not fully 
understood. Their values are therefore best determined 
phenomenologically. Luckily the quantities R{uj) and 
C{uj) depend only on the normal state band-structure 
model, and on the maximum of the mean-field d-wave 
energy gap A. At the current mature state of cuprate ex- 
perimental research, the information we require is there- 
fore available with adequate accuracy. 

Fig. [DJa) plots R{uj) for ui < flo, where ilo is the mini- 
mum two-particle excitation energy, with two very differ- 
ent band structures. Fig. [IJb) plots C{uj). Although the 
association of C with Berry phases is not strictly appro- 
priate away from the adiabatic limit, it is appropriate to 
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FIG. 1: (a) RiiiS) and (b) C{lS) for two difference band struc- 
tures. f2o ~ 70 meV in both models. The property that R(io) 
approaches 1 for a; — + Slo is a universal feature which applies 
for any band structure model. C{lo) is finite and large in both 
models. 



retain the terminology because C(lS) has weak frequency 
dependence up to quite high frequencies. Model I in Fig. 
[T]refers to the band structure model obtained by Norman 
et alw^ by fitting to angle-resolved-photoemssion data^ 
and used by us in our previous work^. Model II refers 
to a model used by Hao and Chubuko\J^ in a related re- 
cent work. The main difference between these two mod- 
els is that the nearest-neighbor hopping t is almost 2.5 
times bigger in Model II than that in Model I. Since the 
maximum energy gaps in both models are fitted to be 
around 35 — 50 meV in the underdoped regime, the ratio 
of the energy gap to the kinetic energy is much smaller 
in Model II than that in Model I. This difference is re- 
sponsible for the generally weaker interchannel coupling 
effects seen in both the R{io) and C(w) in Model II as 
shown in Fig. [T]. We note however that R(lo) ^ 1 as 
approaches JIq in both models, and indeed as we discuss 
below, for any band-structure model. Although Model I 
seems to be much more sophiscated and better fitted with 
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experimentsii, we study both models to ease comparison 
with Ref. [To| and to provide some quantitative indica- 
tion of the sensitivity of our conclusions to very-different 
band-structure model choices. 



III. WEAK INTERACTION LIMIT 

Weak interactions can compete with band energies for 
uj smaller than but close to the minimum two-particle 
excitation energy fig because the energy denominator in 
Eq. [5] becomes small when k approaches p, the point 
at which E{k) -\- E{k') is minimized. By expanding the 
energy denominator around its minimum we find that all 
elements of the quasiparticle response function diverge 
logarithmically as u; ^ fio from below: 



Xll,gp(^ = ^o) 

Xii,qp{^ = ^0 ) 



flip) ln(l - w/f^o) 
flip) I/O ln(l - w/l^o) 
/i(p)/4(p)i^oln(l-c^M). (7) 

Here vq is the joint density of states per site just above the 
minimum excitation energy. Note that all off-diagonal 
elements of the quasiparticle response matrix approach 
the geometric mean of the corresponding diagonal ele- 
ments, so that the response matrrix becomes singular 
and Riuj = O^) ^1. As we explain below, the famil- 
iar BCS-like weak-interaction result for the spin-exciton 
resonance frequency is qualitatively altered as a conse- 
quence. 

When the interchannel coupling is neglected, the spin 
exciton energy ujex is obtained^ by inserting Eq. ([7]) in 
Eq. (O to obtain 
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An undamped resonance appears just below the two- 
particle continuum edge for an infinitesimally small pos- 
itive Vs- To account for the interchannel coupling for lo 
near we rewrite Eq. ([3]) as 

(-X44,gp(w) - Diijj)Vs)i-X\\,qpi^) + OiuS)V^) = 



where 



-D(t^) = (Xll,gp(w)X44,gp(w) - X?4,gp 



M)- (10) 



Note that, although each matrix element diverges loga- 
rithmically the geometric product structure of the matrix 
implies that -D(w) diverges like In, not like In^. After 
some algebras, we obtain 
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(11) 



If > 0, must be greater than a critical value in order 
to have a resonance and ujjnsr > ^ex in agreement with 
Ref. p^ . If < 0, LOiNSR ^ i^ex and an infinitesimal 
positive Vs is again sufficient to guarantee a resonance 
below the two-particle continuum. 



IV. STRONG INTERACTION LIMIT 

Strong interactions drive the system close to the an- 
tiferromagnetic instability, i.e. into the regime which 
appears to be relevant for optimally doped and moder- 
ately underdoped cuprates. Given the instability crite- 
rion, K^iuj = 0) = K^piuj = 0)-Vs with Vs ^U + 2J 
we see that either on-site Coulomb interactions U, or 
Heisenberg near-neighbor spin-dependent effective inter- 
actions J, or a combination of the two, could be respon- 
sible. Since, as we mention again below, experiments 
imply that K'^iuj — 0) is around an order of magnitude 
smaller than K^piuj = 0), the value of Vs is tightly con- 
strained on a relative basis. It has been argued^*' that uj^x 
should match cUres accurately in a relative sense in this 
limit because the frequency which appears in the Berry 
phase coupling term is small. As we will prove in detail 
later, this argument is not necessarily correct because the 
leading frequency dependent term in the low-frequency 
expansion of Eql3] varies as w^, and all terms, includ- 
ing the Berry phase coupling term, should be retained. 

Let's consider the case in which the system is close 
to the antiferromagnetic instability so that the INSR 
frequency is small. We expand the quasiparticle- 
response dependent quantities at low frequencies as fol- 
lows: Xll,9p('^) ~ -(-Rg + R2^^), Xii^qpi^) ~ -(^0 + 

i?2W^), i?(w) « i?i4u;^, and Cicv) ~ C. Detailed expres- 
sions for the positive-definite coefficients Rq'^ , i?2"^ and 
i?i4 can be found in Ref. 0. Ri4 and C are related 
by C = Rii/V where V is the low frequency limit of 
Xi4,gp(^)/'^- The spin-exciton energy can be obtained 
by setting T4xii,9p(^) = ^1- We find that 
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where S is defined by the second equality. The above 
expression is valid close to the antiferromagnetic insta- 
bility, i.e. for 6 =—f 0"*". When the interchannel coupling 
is included the corresponding resonance energy luinsr 
expression is obtained by solving Eq. [3] retaining terms 
up to order w^. After some algebras we find: 

i^INSR 1 _ ,,„^ 

7 (13) 



where 



and 



OJe 



VI + 71 + 72 



71 



(l-<5)(l + yX) 



- 1 



72 = 



-5iRt - RtRu)V^ 



(14) 



(15) 



(1 + V^Rl)RiVs 

Since 72 oc 5 — > near the spin-density instability, the 
discrepancy between loin sr. and ujex is mainly associ- 
ated with the parameter 71. We comment further on 
this point in the next section. Clearly the value of 71 
depends strongly on model parameters. 
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V. DISCUSSION 

As emphasized in our previous worki, because we 
have three interaction parameters {U, V, J) but only two 
experimentally-determined properties, the anti-nodal 
gap and the INSR frequency, experiment does not 
uniquely determine the effective interaction model. The 
two combination of interaction parameter combinations 
that are tightly constrained by experiment are the spin- 
interaction Vs = U + 2J (INSR experiment) and the pair- 
potential Vpair = 3J/2 — 2V (ARPES and other measures 
of the d-wave antinodal gap). We have argued previously 
that the V = 0, J ( J-scenario) possibility (pairing 
driven by effective Heisenberg interactions) is more likely 
than that V ^ 0, J — (T^-scenario) possibility (pairing 
driven by attractive scalar near- neighbor interactions). 
Our argument appealed to the Occam's razor assump- 
tion that one type of effective interaction is likely to play 
a dominant role in the low-energy physics. We then ob- 
served that the J-scenario can account for both Vg and 
Vpair scales whereas the ^-scenario requires in addition a 
finely-tuned on-site effective interaction. ([/ would have 
to be at least 350meV in model I and SOOmeV in model 
II if the ^-scenario was adopted to fit the experimental 
data^»i2i.) It is in fact clear from the outset that the near- 
neighbor scalar interaction V cannot succeed in account- 
ing for both scales because it contributes to Vpair but 
not to Vs = U + 2 J. On the other hand, the J-scenario 
is not trivially guaranteed to explain both experimental 
energy scales. This fact does appear to us to be strongly 
suggestive. 

Eqs. (|13I14I15[) demonstrate that interchannel cou- 
pling can alter the resonance frequency even in the limit 
Kg and that the resonance frequency always de- 
pends on both V^ and V^^ . In the J-scenario we favor is 
small and positive, and 71 tends to be negative in models 
which are close to the antiferromagnetic instability. For 
this reason 7 is slightly larger than 1 in both models we 
have examined. For Model I, we find 7 « 1.1 while for 
Model II we find 7 w 1.04. The approximate coincidence 
between LUex and loinsr for some model parameters does 
not imply that particle-hole and particle-particle modes 
are weakly coupled in the elementary excitations. In fact, 
it is evident in Fig. ([1]) that C{uj) is weakly frequency 
dependent and remains comparable to 1 in both models, 
which already guarantees the strong interchannel cou- 
pling as discussed above. The close match between uj^x 
and wiNSR only reflects that the choice of interaction pa- 
rameters leads to not only K^{uj) (i.e. S ^ 0) but 
also the large ratio of K^{uj)/ Ks{u!) at low frequency. For 
realistic experimental situations in which the resonance 
frequency is always a substantial fraction of the antin- 
odal gap, this ratio is always larger than approximately 
3, which does favor ujinsr ~ ^ex- This analysis proves 
that the ratio of ujiNSR./^ex can not serve as a reliable 
criterion to judge the importance of the interchannel cou- 
pling. We have referred to this resonance as a magnetic 
plasmon to suggest this aspect of its character, namely 



that it represents the coupled quantum fluctuations of 
canonically conjugate degrees of freedom as long as |C| 
is comparable to 1. Our conclusions on this point arc in 
disagreement with those of Ref. [l3| • 

In Ref. 3 we have argued that the quantum zero- 
point energy of cuprates makes a negative contribution 
to the superfluid density of cuprates, and that it there- 
fore contributes to the decrease of the critical tempera- 
ture in the underdoped regime. The basic idea is that 
superconductivity is weakened by phase gradients, al- 
lowing the antiferromagnetic instability to creep closer 
and reducing the quantum zero-point energy associated 
with spin-fluctuations. Similar physics occurs in bilayer 
quantum Hall superfluids^^'^'^. Microscopic caclulations 
suggest that in cuprates this effect is large enough to ac- 
count for the decreasing superfluid density, in contrast 
to the situation in typical superconductors in which cor- 
relation effects are largely independent of the configura- 
tion of the superfluid condensate. It will be interesting 
to examine how this idea plays out when applied to the 
recently-discovered iron pnictidc superconductors. Cur- 
rent experiments and theories seem to favor extended-s- 
wave gap symmetry. An INSR has been observe d^"^'^^i^^ 
in the pnictide materials with a phenomenology which 
is surprisingly similar to that of the cuprates. Because 
the the ratio 2A/fcBTc is found to be remain close to 
the BCS predcitioi*ii, it seems that the supression of su- 
perfluid density by correlations is less important in this 
material. To account for this difference within our theory, 
we will need to understand why phase gradients in the 
superconducting condensate have a weaker effect on spin- 
fluctuations and on the correlation energy. The multi- 
band character and the peculiar Fermi surface pockets of 
pnictides would seem likely to produce this result. The 
large effect in cuprates rests strongly on the importance 
of fluctuations near the wavevector which connects anti- 
nodal momenta for spin- response functions, a coincidence 
that is not likely to be repeated in iron pnictide supercon- 
ductors. A recent RG analysisi^ which showed that the 
particle-hole and particle-particle channels are decoupled 
at energy scalies below the Fermi energy Ep supports 
this view. A more detail calculations for the iron pnic- 
tides using an appropriate multiband model which aims 
to resolve this issue is currently in progress. 

Based on these analyses, we conclude that Berry phase 
coupling between the particle-hole and particle-particle 
channels can never be ignored in any RPA-type theory 
of the inelastic neutron scattering resonance mode for 
cuprates. As a result, the inelastic neutron scattering 
resonance mode observed in cuprates should be bettter 
interpreted as a magnetic plasmon than a spin exciton. 
The importance of coupling in the pnictides has not yet 
been adequately analyzed. Nevertheless, the spin res- 
onance mode observed in the inelastic neutron scatter- 
ing measurements could be a purer spin exciton, and the 
magnetic and superconducting instablities less coupled in 
the pnictide case. 
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